Thermal rectification effect of an interacting quantum dot 
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We investigate nonlinear thermal transport properties of a single interacting quantum dot with 
two energy levels tunnel-coupled to two electrodes using nonequilibrium Green function method 
and Hartree-Fock decoupling approximation. In the case of asymmetric tunnel-couplings to two 
electrodes, for example, when the upper level of the quantum dot is open for transport, whereas the 
lower level is blocked, our calculations predict a strong asymmetry for the heat (energy) current, 
which shows that, the quantum dot system may act as a thermal rectifier in this specific situation. 

PACS numbers: 74.25.Fy,72.10.Bg,73.63.Kv,73.23.Hk 
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Thermal transport properties of various nanoscale de- 
vices, e.g., a semiconductor quantum dot (QD) and a 
molecular junction, have recently gained considerable 
attention^ Most theoretical work on thermoelectric effect 
of a QD attached to external leads have so far focused 
on the linear thermopower in sequential^ and cotunnel- 
ing regimes,— and even in the strongly correlated Kondo- 
type transport regime4^ In addition, quantum interfer- 
ence effect in the thermopower has also been studies in a 
QD-ring geometry£>£ A great amount of investigations of 
heat conduction through molecular chains connecting two 
electrodes at different temperature have been devoted to 
building a thermal rectifier, an analog of charge current 
rectifier first proposed by Aviram and Ratner in molec- 
ular transistor^ which transport heat current (phononic 
current there) efficiently in one direction of the applied 
temperature difference between two electrodes, but block 
it in the reverse direction^ The main ideal to raise a 
charge or thermal rectifier there is to make asymmetrical 
couplings between the device and electrodes, which is vi- 
able due to advances in the fabrication of nanostructures. 

A very recent experiment by R. Scheibner et.al. has 
reported that a semiconductor QD with two energy lev- 
els can also behave as a thermal rectifier under some 
conditions^ To achieve asymmetric tunnel-couplings to 
two electrodes, they tuned high orbital momentum states 
of the QD as the active transport channel and applied 
high in-plane magnetic fields to quench the channel. This 
experimental results are quite interacting and very useful 
in thermoelectric application since the predicted rectifier 
is performed only in the electronic system and without 
the need for coupling to the phonon system. Stimulated 
by this experiment, we theoretically study, in this pa- 
per, the nonlinear heat transport through an interacting 
QD with two energy level asymmetrically coupling to two 
electrodes when a temperature difference is applied be- 
tween the two electrodes, and examine in detail the rea- 
son and condition of appearance of thermal rectification 
effect in this system. 

The rest part of this paper is as follows. In Sec. ??, we 



introduce the model of a single QD with two energy lev- 
els and give the main theoretical formulation employed 
in this paper. In Sec. ??, we carry out numerical calcula- 
tions and some detailed analyses with different coupling 
cases. Finally, a brief summary is given in Sec. ??. 

In this paper, we consider a single interacting QD with 
two energy levels tunnel-coupling to two normal leads, 
as illustrated in Fig. [TJ The model Hamiltonian of this 
system can be written as: 



H = H L + H h 



Hr 



H, 



(1) 



where H v (rj = L,R) describes noninteracting electron 
baths in the left and right leads, respectively: 



Hv ~ E £, ) kC »7k C r/k' 



(2) 



with c^ k (c,jk) being the creation (annihilation) operator 
of an electron in lead rj with momentum k, and energy 
measured from the Fermi energy of the electrode. The 
chemical potentials of the left and the right leads ^l/r 
are associated with the bias voltage eAV — fiL ~ MR) 
which is applied symmetrically between the two elec- 
trodes, i.e., fiL = —H>r = eAV/2. Similarly, we as- 
sume that a temperature difference AT is also applied 
symmetrically between two leads, Tl = T + AT/2 and 
T R = T - AT/2 (T is the initial temperature of both 
electrodes). In equilibrium, we have AT = AV = 0. Hjj 
describes the QD with two energy levels: 



Hr 



E 

J=l,2 



^ C 1 C 1 C 2 C 2> 



(3) 



where 



,•) is the annihilation operator for a spinless 



electron in the jth level (j = 1,2), with Sj being the cor- 
responding level energy. The second term represents the 
Coulomb interaction U between two energy levels. The 
final term in Eq. (JTJ), Hj, describes the tunnel-coupling 
between the levels and the electrodes: 



Hr = 



H.c. 



(4) 
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Here, we ignore spin degree of freedom for simplicity. 
Throughout, we will use units with h = kg = e = 1. 




FIG. 1: Schematic diagram for a single QD with two energy 
levels coupling to two electrodes. 



We consider a nonlinear system with a finite temper- 
ature difference AT and a bias voltage AV applied be- 
tween the two leads, which will arouse electronic current 
J c and heat current J h flowing through the QD. Apply- 
ing the Keldysh's nonequilibrium Green function formal- 
ism, one can obtain: 10 



r = — I dui T(ui) [jxM - f R (u)], 



(5) 



J h = ^ J dw(w - AV) T{u) [f L (cj) ~ f R (uj)}, (6) 

where fr,(tj) = (e' w_M "" T " + 1) _1 is the Fermi distri- 
bution of lead rj. The transmission coefficient, T(u>), 
can be expressed using the Fourier transformation of 
the retarded and advanced Green functions of the QD, 
[G^(t,t')] jf = TiO(±tTt')({cAt) 7 cUt>)}), as: 



T(uj) = Tr 



TrT 



Li- R 



[CT(w) - G»] 



(7) 



with the constant level- width function [T]^/ = 
2TTQ v V V jV V ji {g v is the density of states of lead rj). So T v jj 
is the bandwidth of the jth level of the QD due to the 
coupling with lead rj. The retarded Green function can be 
obtained using the standard equation of motion method 
with help of a decoupling approximation. In this paper, 
a Hartree-Fock decoupling scheme is employed and is suf- 
ficient for the present purpose since we take no account 
of the Kondo effect and assume a higher temperature 
in our calculations, T 3> Tk (Tk is the Kondo temper- 
ature). For example, we derive the diagonal retarded 
Green function as: 



- e 3 - U(l - nj) 



W = 

H(w) = (uj - e 3 )(uj - £j -U) 
-E^- £j --t/(l 



nj)], 



(8) 
(9) 



with j = 1(2) if j = 2(1). = -(i/2)(T Lj + T Rj ) is 
the retarded self energy of the QD due to coupling to 



the electrodes, and rij = (cjcj) is the occupation num- 
ber of electrons in the jth level, which can be calcu- 
lated self-consistently through the lesser Green function, 
[G<(t,f)] jjt =i(c],(t%(t)), as 



(10) 



The lesser Green function is related to the retarded Green 
function by the Keldysh function G< = G r H<G a and 
S j> = ^ljMuj) +r Rj f R (u)]5 jj >. 

Based on these equations (f5|)- ([T0l) . we can study the 
thermoelectric properties of the QD in the two-terminal 
configuration, i.e., the heat transport through the QD 
when a temperature difference AT is applied between 
the two electrodes but the charge transport is vanishing 
due to an open circuit condition. This external applied 
finite temperature difference will cause diffusion of elec- 
trons cross the device from one terminal to the other and 
generate accumulation of electrons at one of the termi- 
nals under the condition of open circuit. Such accumu- 
lation gradually builds a nonzero bias-voltage across the 
QD to counteract the diffusion of electrons until a steady 
state is finally reached. Numerically, this temperature- 
difference-induced bias voltage must be obtained by solv- 
ing the equation 



J c = 0. 



(11) 



The current formula, Eq. ([5]), is relevant to the electron 
occupation numbers, n\ and ri2, through the Green func- 
tions, Eqs. ijHJ) and Meanwhile, the electron occupa- 
tion numbers are reversely dependent on the temperature 
difference and induced bias voltage through Eq. pTJj) . 
Therefore, for a given temperature difference, AT, we 
have to solve both Eqs. (|10p and (TTTj) simultaneously to 
obtain a self-consistent result for the induced bias volt- 
age, AV, and then use it in Eq. © to calculate the heat 
current through the QD. 

Based on the scheme indicated above, we investigate 
the nonlinear thermoelectric transport through an in- 
teracting QD with two energy levels to search for the 
possible thermal rectification mechanism. As motivated 
by the recent experimental results^ we focus our stud- 
ies on the two cases of system: the symmetric tunnel- 
coupling case, T R j = r^ , and the asymmetric case with 
a blocked lower level, T^i = 0. In the following numeri- 
cal calculations, we choose Tli = T l 2 =T as the energy 
unit. We also set the temperature of equilibrium situa- 
tion as T = 4.0 r, the inter-level Coulomb interaction as 
U = 100 r, and fix the energy gap between the two levels 
as E2 — s\ = 10 T. 

We show our main results in Fig. [21 the calculated 
nonlinear heat current J h as a function of temperature 
difference AT with different couplings between the two 
levels and electrodes. It is natural that for the symmetric 
coupling system (1), the heat current exhibits symmetric 
behavior, i.e., monotonously increases with increasing the 
temperature difference in both directions, AT > and 
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AT < 0. In contrast, a strong asymmetric behavior is 
found for the asymmetric systems (2) and (3), whose heat 
current increase with the temperature difference only in 
the positive direction (AT > direction), but is blocked 
in the opposite direction (AT < 0), indicating the occur- 
rence of a thermal rectification effect. It seems that this 
thermal rectification effect is a direct result of the asym- 
metric tunnel-coupling arrangement of the QD. However, 
the situation is not so simple because our calculations ex- 
hibit no thermal rectification effect for the system with 
blocking the higher state instead of the lower state [see 
Fig- SJa) below]. Recently, an experimental observation 
of thermal rectification effect of a semiconductor QD has 
been reported by R. Scheibner, et.al.M- In their measure- 
ment, they set up some high orbital momentum states of 
the QD to play an important role in heat transport and 
applied high in-plane magnetic field to control the cou- 
pling strength between these states and the electrodes. 
They measured the gate- voltage-dependent thermopower 
of the QD and found that it changes to an asymmetric 
line shape when the lower state is blocked by applying a 
high magnetic field, which reveals the thermal rectifica- 
tion effect. Our calculation results are in agreement with 
this experimental inferences. 
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FIG. 2: Heat Current vs. temperature difference with the 
system parameters ei = —6, £2 = 4, for symmetric coupling 
case (1) r_Ri = 1, r_R2 = 1, asymmetric cases (2) Tri = 0, 
Fr2 = 1, and (3) r^i = 0, Tr2 = 10. Inset: rescaling plot 
for the asymmetric case (2). The other parameters are T = 4 
and U = 100. All are in unit of T. 

In order to interpret this effect, we plot the self- 
consistently calculated electron occupation numbers, ri\ 
and U2-, in the two levels as functions of the applied tem- 
perature difference in Fig. [3] It is visible that 71.1(^2) of 
the symmetric system always increases (decreases) with 
increasing the temperature-difference in both directions, 
and the total electron number n\ + n-i also increases. 
It means that with increasing temperature difference (in 
both directions), more and more electrons participate 
in the energy transferring processes, thus resulting in a 
monotonous rise of the heat current. If the lower state 



is blocked, we find that to fulfill the open circuit condi- 
tion, the induced bias-voltage AV has a negative slope 
with respect to AT for the specific system parameters 
(not shown here). Therefore, in the region of AT > 0, 
the negative AV indicates a lowering chemical potential 
of the left lead /i£ with a heightening temperature Tl, 
implying a lesser occupation possibility /l(ei) for elec- 
trons in the left lead to enter the QD, finally resulting 
in a lesser n\ and corresponding a larger ri2 due to the 
strong Coulomb correlation. In contrast, in the region 
of AT < 0, a larger Fermi distribution makes n\ 

increase more quickly than that of the symmetric sys- 
tem with AT arising, and 712 decrease nearly towards 0. 
Since the upper level (level-2) is the only open transmit- 
ting level in the asymmetric system, it is natural that the 
large 712 in the positive direction of AT leads to a large 
heat current, while the vanishing ri2 in the negative di- 
rection is responsible for a suppressed electron transfer, 
i.e., indicating the occurrence of a thermal rectification. 
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FIG. 3: Electron occupation numbers vs. temperature differ- 
ence with the same system parameters as in Fig. [2] 

However, the situation is different for the upper state 
blocked case (Tr2 = 0) as shown in Fig. 2J in which 
the calculated heat current J h , and electron occupation 
number of the lower level, n\ , are plotted as functions of 
temperature difference. It is found that the heat current 
increases in both directions, and there occur no thermal 
rectification effect. This is physically because the lower 
electronic state is now the open transmitting level, thus 
the occupation number of this level, m, determines the 
heat current. From Fig. [31(b) , it is clearly to see that 
n\ keeps a relatively large value at the whole region of 
AT. That is to say that there are sufficient numbers of 
electrons to carry energy transferring in both directions. 

In above discussions, the chemical potentials of both 
electrodes are between two energy levels at equilibrium. 
One may think that it is the special arrangement of the 
level energies which decides the thermal rectification ef- 
fect. In fact, to raise a thermal rectification effect, one 
needs to suppress the occupation number of transmitting 
level in one direction of AT. For a two-level system, the 
lower level always holds a large occupation number. Thus 
if and only if the lower level is blocked, one can observe 
the thermal rectification. Besides, a large Coulomb inter- 
action U is also necessary for the rectification, because it 
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FIG. 4: Heat Current J h , electron occupation number of the 
lower level n\ vs. temperature difference for symmetric cou- 
pling case (1) r_Ri = 1, F_R2 = 1, upper state blocked case 
(2) F m = 1, F R2 = 0, and (3) T m = 10, r R2 = 0. Other 
parameters are the same as in Fig. 2. 

leads to large repulsion between n\ and n,2, and thus can 
greatly suppress the occupation number of the upper 
level. To verify the these conclusions, we make calcula- 
tions for two cases where two energy levels of the QD are 
both above or below the chemical potentials of the two 
electrodes at equilibrium, as shown in Fig.[5l It gives the 
expected result that a thermal rectification effect occurs 
for the blocking lower level cases. The rectification effect 
m Ref.0 indeed corresponds to the case where two levels 
are both above the chemical potentials with blocking the 
lower state. We also calculate the nonlinear heat trans- 
port without the inter-level Coulomb interaction and find 
that there occurs no rectification for all symmetric and 
asymmetric systems (not shown here). 



FIG. 5: Heat Current J h vs. temperature difference for cases 
with two energy levels both above (a) or below (b) the chem- 
ical potentials of two electrodes, with parameters (a) ei = 4, 
£2 = 14, and (b) ei = —16, £2 = —6. Lines (1),(2),(3) are re- 
spectively corresponding to the symmetric coupling case (1) 
r_Ri = 1, r_H2 = 1, lower state blocked case (2) Tri = 0, 
r_H2 = 1, and higher state blocked case (3) F_ri = 1, Fr2 = 0. 

In summary, we investigated the nonlinear thermoelec- 
tric effect of a two-level quantum dot system with inter- 
level Coulomb interaction. We calculated the nonlinear 
heat current through the quantum dot when a temper- 
ature difference is applied between the two electrodes, 
and found that in some cases the heat current displays 
a strong asymmetrical behavior, which indicates a ther- 
mal rectification effect. We analyzed further this effect 
in detail, and ascribed it to the asymmetrical couplings 
between two electrodes and two levels: only when the up- 
per level is open for transport and meanwhile the lower 
level is blocked, the rectification effect may occur. Be- 
sides, our numerical discussions shown that the strong 
Coulomb interaction is necessary for the occurrence of a 
rectification effect. 
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